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ABSTRACT 


The  ability  to  predict  torsional  instability  in  the 
early  stages  of  design  can  have  important  consequences 
on  the  design  of  both  conventional  and  high  performance 
ships.  This  thesis  develops  fast  approximate  methods  of 
torsional  buckling  analysis  for  particular  application 
in  the  concept,  feasibility,  and  preliminary  stages  of 
ship  design. 

Two  simplified  models  of  stiffeners  commonly  used  in 
ship  construction  were  presented.  The  first  was  an 
ideal  I  section  st if fener.  The  second,  a  flat  bar 
stiffener,  which  was  included  in  this  analysis  because 
of  its  favorable  product ional  properties  and  wide 
useage.  These  models,  with  initial  imperfections,  were 
subjected  to  axial  compressive  end  loading  and  the 
resultant  behavior  analyzed.  The  approach  of  the 
analysis  of  this  thesis  was  an  application  of  the 
energy  method  to  determine  the  critical  buckling  stress 
instead  of  the  more  commonly  used  equilibrium  approach. 
Both  beam  theory  and  thin  plate  theory  were  used  in 
energy-work  relationships  in  these  derivations.  Then 
the  first  yield  load  was  determined  utilizing  the 
Perry-Robertson  approach.  Integrated  into  all  phases  of 
this  exploration  was  the  concept  of  the  geometric 
imperfections  of  the  stiffeners. 
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NOMENCLATURE 

A   Nondimensional  parameter 

a   Length  of  stiffener  between  transverse  supports 

b   Uniform  st if fener  spacing 

b   Plating  effective  width 

e 

B   Angle  of  rotation  of  the  st if fener  about  the  toe 

* 
B    Initial  imperfection  of  the  st if fener  in  terms  of  a 
rotation  angle  about  the  toe 


B  Amplitude  of  the  sinusoidal  function  for  B 
o 

*  * 

B  Amplitude  of  the  sinusoidal  function  for  B 
o 

C  Rotational  spring  constant  (moment/length)  of   the 
supporting  plating 

C  Longitudinal  warping  constant  of  stiffener 
w 

D  Flexural  rigidity  of  plating 

D  Flexural  rigidity  of  web  plating  of  stiffener 
w 

d  Depth  of  flat  bar  stiffener;   overall  depth  of  Tee 
st  if fener 

d  Depth  of  stiffener  to  midthickness  of  flange 

c 

d  Depth  of  web 
w 

E  Young's  Modulus  of  the  material 


e"    Membrane  strain 
x 

f    Width  of  stiffener  flange 
w 

G   Shear  Modulus  of  the  material 

h   Height  of  neutral  axis  of  plate-st if fener 
combination  from  midplane  of  plating 

I    Effective  vertical  moment  of  inertia  of   stiffener 
and  associated  effective  width  of  plating 

I    Polar  moment  of  inertia  of  stiffener  about  toe 
P 

I    Vertical  moment  of  inertia  of  stiffener  alone  about 
t 

toe 

I    Moment  of  inertia  of  stiffener  about  web  plane 

2 

J   St.  Venant's  torsion  constant  for  stiffener 

L   Reduced  slenderness  ratio 

M(x)    Moment  developed  in  the  stiffener  by  the  total 
torque 

m   Mode  number 

P   flxial  end  load 

s   Height  of  stiffener  shear  center  above  toe 

£     Elastic   axial   torsional   buckling    stress    in 
CR 

stiffener  (critical  stress) 


£    Elastic   axial  torsional  buckling  stress  limit   of 
CL 

stiffener  (classical  stress) 


5"   Axial  stress  in  stiffener 

e 


£   Axial  stress  component 
x 

(%"  )       Maximum  compressive  stress  developed   by 
x  MflXT 

the  torque 


s*   Tensile  yield  stress  of  material 

y 


T   Amplitude  of  the  sinsoidal  torque  function 
o 

T   Total  torque  developed  in  the  stiffener   structure 

by  the  loading 

T   St.  Venant's  torque  developed  in  the  stiffener 
SV 


T   Warping  torque  developed  in  the  stiffener 
w 


t  Plate  thickness 

t  Stiffener  flange  thickness 
f 

t  Stiffener  web  thickness 
w 

u  Displacement  in  the  longitudinal  (x)  direction 

V  Total  strain  energy  of  the  structure 

v  Sideways  flexure  (y)  of  stiffener 

v  Initial  imperfection  of  stiffener  in  horizontal  (y) 
direction 

W  Total  work  of  the  applied  force 

w  Vertical  flexure  (z)  of  stiffener 


* 
w    Initial   imperfection   of   the   stiffener   in   the 
vertical  (z)  direction 

z   Height  of  stiffener  centroid  above  toe 

"V   Poisson* s  Ratio 


Note:   x,y  or  z  subscripts  indicate  partial  derivatives 

with   respect   to  those  coordinates  .   These   are   used 

*      *      * 
particularly  with  v,  v  ,  w,  w  ,  B,  B  . 
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INTRODUCTION 

The  subject  of  this  thesis  is  to  develop  a 
simplified  analytical  method  of  predicting  torsional 
buckling  of  stiffeners  used  in  ship  structures. 
Torsional  buckling  or  instability  is  the  phenomena  by 
which  a  column,  at  a  certain  axial  load,  fails  suddenly 
in  a  combined  mode  of  twist  and  lateral  bending  of  the 
cross  section  (reference  7).  Figure  (0-la)  (from 
reference  6)  illustrates  the  failure  mode  referred  to 
here  as  torsional  buckling  -  in  this  case  for  a  flat 
bar  stiffener.  Figure  (0-lb)  illustrates  the  failure 
mode  of  the  tripping  phenomenon  which  is  not  covered  in 
this  analysis.  Figure  (0-2)  (reference  6)  is  the 
experimental  load  deflection  curve  for  the  failure  of  a 
flat  bar  stiffener  under  an  axial  load.  Point  ft  of  the 
solid  (experimental)  curve  of  figure  (0-2)  is  the 
ultimate  strength  load.  As  can  be  seen  from  figure  (0- 
2),  the  curve  proceeds  in  a  downward  fashion  -  the  bar 
looses  stiffness  and  fails  from  this  point. 

The  ability  to  predict  torsional  instability 
in  the  early  stages  of  design  can  have  important 
consequences  on  the  design  of  both  conventional  and 
high  performance  ships.  It  is  the  intention  of  this 
thesis  to  develop  fast  approximate  methods  of  torsional 
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Figure  (0-la)   Torsional  Buckling 


Figure  (0—1 b)   Tripping 


Figure  (0-1)  Characterization  of  Torsional  Buckling 
and  Tripping  in  a  Flat  Bar  Stiffener 
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Figure  (0-2)  Experimental  Load  Deflection  Curve 
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buckling  analysis  for  particular  application  in  the 
concept,  feasibility,  and  preliminary  stages  of  ship 
design. 

To  this  end,  two  simplified  models  of  stiffeners 
commonly  used  in  ship  construction  will  be  presented. 
The  first  is  an  ideal  I  section  stiffener.  The  second, 
a  flat  bar  stiffener,  which  is  included  in  this 
analysis  because  of  its  favorable  product ional 
properties  and  wide  useage.  These  models,  with  initial 
imperfections,  will  be  subjected  to  axial  compressive 
end  loading  and  the  resultant  behavior  analyzed.  The 
approach  of  the  analysis  of  this  thesis  will  be  an 
application  of  energy  methods  to  determine  the  critical 
buckling  stress  instead  of  the  more  commonly  used 
equilibrium  approach.  Both  beam  theory  and  thin  plate 
theory  are  used  in  energy-work  relationships  in  these 
derivations.  Then  the  first  yield  load  is  determined 
utilizing  the  Perry-Robertson  approach  as  described  in 
references  (9,14).  Integrated  into  all  phases  of  this 
exploration  will  be  the  concept  of  the  geometric 
imperfections  of  the  stiffeners. 

It  is  the  intention  of  this  thesis,  using  the  above 
mentioned  approach,  to  lay  groundwork  which  would  have 
the  potential  for  future  expansion  to  include  mode 
interaction  -  to   be  able  to  easily   analyze   combined 
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loadings  etc. . 

Particular  emphasis  is  placed  here  upon  inclusion 
of  the  geometrical  initial  deformities  of  the 
stiffeners  which  occur  during  the  manufacture  and 
through  day  to  day  usage  of  the  ship.  It  is  well  known 
that  these  geometric  imperfections  can  drastically 
reduce  the  strength  of  the  structural  member.  But 
current  design  standards  (references  10,11  etc)  have 
not  tried  to  quantify  these  values.  Some  design  codes 
(example:  reference  10)  have  specified  tolerence  limits 
of  various  geometric  parameters  which  must  be  met 
before  the  design  formulae  are  considered  valid. 
Granted,  in  the  concept,  feasibility  and  preliminary 
stages  of  ship  design  where  the  ship  only  exists  on 
paper,  the  designer  would  have  no  knowledge  of  the  end 
resulting  imperfections  of  the  actual  ship.  Hopefully, 
in  the  future,  with  the  advent  of  better  production 
control  in  shipyards  and  the  incursion  of  the  computer 
into  ship  production  facilities,  this  information  will 
be  collected  and  statistically  analyzed. 

In  addition,  once  the  ship  is  manufactured,  these 
formulae  will  provide  a  quick  and  easy  evaluation  of 
the  actual  critical  buckling  stress  in  terms  of 
physical  parameters.  Thus,  this  analysis  has  made  a 
special   effort   to  define  the  imperfections  in   easily 
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measured  terms.  This  analysis  measures  the  horizontal 
or  vertical  movement  of  the  shear  center  of  the 
stiffener  at  its  point  of  greatest  deflection  and 
assumes  a  simple  sinusoidal  distribution.  Knowledge  of 
the  web  depth  easily  converts  the  imperfection  into  a 
rotational  angle.  Here,  the  stiffener  is  considered  to 
rotate  about  its  toe  or  base.  In  other  words,  the  base 
plating  is  considered  much  larger  and  stiffer.  Figure 
(0-3)  illustrates  this  simplified  concept  for  a  Tee 
st  if fener. 


Figure  (0-3)  Tee  Stiffener 

Currently  the  method  of  handling  the  stiffener 
failure  problem  is  incorporated  into  the  use  of  factors 
of  safety.  But,  with  recent  attempts  to  further 
optimize   ship  structures  -  in  particular   to   minimize 


IB 


the  amount  of  weight  of  a  ship  that  must  be  devoted  to 
its  structure  -  better,  yet  reliable,  methods  must  be 
devised.  Remember,  with  a  ship,  less  weight  in  the 
structure  equates  to  more  weight  for  the  pay load  be  it 
weapons  for  a  military  ship  or  cargo  for  a  civilian 
ship. 

This  thesis  concludes  with  a  comparison  between  the 
formulae  derived  from  the  simplified  models  presented 
here  and  published  formulae,  published  finite  element 
analysis,  and  experimentally  derived  results. 
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1.  DESCRIPTION  OF  MODELS 

Torsional  buckling  or  instability  is  characterised 
by  a  twisting  of  the  stiffener  about  its  line  of 
attachment  to  the  plating.  This  deformation  pattern 
involves  both  sideways  and  vertical  flexure  (v, w)  and 
rotation  (B)  of  the  stiffener  as  shown  in  figure  (1)  as 
described  in  reference  (3). 


y< 


Figure  (1-la)   Coordinate  System 


± 


Y<r 


Figure  (1-lb)  Undeformed        Figure  (1-lc)  Deformed 
Figure  (1-1)   Characterization  of  Torsional  Buckling 
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1.  1  Model  I._ 

The  first  model  considered  is  a  simplification  of 
an  ideal  I-section.  It  consists  of  an  I-beam  with  zero 
torsional  rigidity.  The  web  of  the  I-beam  transmits 
only  shear  forces.  In  effect,  it  looks  like  two  unequal 
flanges: 


UOOEZZZ    >  >  n 


I-beam  Stress  Distribution 

Figure  (1-2)   Model  I 

This   model   has  been  considered  in  reference   (4). 
However-  reference  (4)  used  an  equilibrium  approach  and 
initial     imperfections    were    considered    as     an 
eccentricity  of  the  axial  load.  The  resulting  equations 
of  the  critical  stress  were  quite  complicated. 

Figure   (1-3)  illustrates  the  geometric  parameters 
of  this  model. 
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Figure  (1-3)  -  Geometric  Torsional  Buckling 

Parameters  for  Model  I 


I   -  Moment  of  inertia  about  the  web  plane  (top   flange 
z 

only) 


I   =  (t  f  )/12 
z     f  w 


w 


s   -  Height   of   shear  center  above  toe   (origin)   (top 
flange  only) 


s  ■  d 

c 

-  Longitudinal  warping  constant  about  shear   center 

of  stiffener  alone  (top  flange  only) 

c  -  a 

w 

-  Vertical  moment  of  inertia  about  toe  (top   flange 

only) 

2  2 

I   »  f  t  (d   +  (l/12)t  ) 
t     w  f   c  f 
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I   -  Polar  moment  of  inertia  about  toe 
P 

1=1+1 

P      t      2 

z  -  Height  of  centroid  above  toe 

2  ■  d 


J  -  St.  Venant's  torsion  constant  (top  flange  only) 
3 
J  =  (f  t  )/3 
w  f 


21 


1.2  Model  II;. 

The   second  model  considered  is  the  simple  flat  bar 
stiffener  illustrated  in  figure  (1-4). 


Figure  (1-4)  Model  II  -  Flat  bar  stiffener. 

Figure   (1-5)  illustrates  the  geometric   parameters 
of  this  model. 
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Figure  (1-5)  -  Geometric  Torsional  Buckling 

Parameters  for  Model  II 


I   -  Moment  of  inertia  about  the  web  plane 

Z 

(web  only) 
3 
I   =  (d  t  )/12 

Z       WW 

s  -  Height  of  shear  center  above  toe  (origin) 
(web  only) 

s  =  d  /2 
w 

C   -  Longitudinal  warping  constant  about  shear 

w 

center  of  stiffener  alone 

C   =  0 
w 

I   -  Vertical  moment  of  inertia  about  toe  (wsb  only) 
t  3 

I   -  (t  d  )/3 
t     ww 

I   -  Polar  moment  of  inertia  about  toe 

P 

I   -  I   +  I 
p    t     z 
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z  -  Height  of  centroid  above  toe 

2  =  d  /2 
w 

J  -  St.  Venant's  torsion  constant 
3 
J  -  (d  t  >/3 
w  w 

D   -  Flexural  rigidity  of  the  web  plate 
w 

3 

Et 
w 

D   = 

w  2 

12(1-  V  ) 
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2.  MODEL  I  -  TORSIONAL  BUCKLING  UNDER  AXIAL  LOADS 

2.  1  Devel_2Brnent  of  the  Strain  Energy  Eguat  i,gn._ 

Since  Model  I  basically  consists  of  two  unequal 
flanges  with  the  top  flange  rotating  about  the  bottom 
flange,  the  torsional  rigidity  is  considered  only  for 
the  top  flange.  Thus  the  strain  energy  equation 
consists  of  three  parts  (reference  1).  The  first  term 
represents  sideways  bending.  The  second  term  represents 
longitudinal  warping.  The  third  term  represents  the 
torsional  rigidity  of  the  top  flange  alone.  The  strain 
energy  for  a  length  a  is: 


r 


(EI  v    +  EC  B    +  GJB  )dx 

0      2  XX        W  XX         X 


Also   from  geometrical  properties  we  have   (see   figure 
1-3)  : 

v  a*  sB 
Thus: 


•a      2  2  2      2 

(EI  s  B    +  EC  B     +  GJB  )  dx 

0      2     XX        W  XX  X 


2  f 

Si* 


a     2       2       2 

|   Ed  s   +  C  )B    +  GJB  1    dx 

0      2         W    XX         X 
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The   following   choice  of  B  is  made,   in  order   to 

meet   certain  boundary  conditions  (ie  B    =0   at   x=0, 

xx 
x=a  )  and  to  keep  B  simple. 


B  ■  B   sin  (mTTx/a) 
o 

Thus  for  Model  I: 

2       2       2  2 

V  ■  1/4  a  B   <m7?Va>  C  E(I  s   +  C  Xm^/a)   +  GJ  1 

O  Z        W 
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2.  2  Sign  Convent  ion 

Before  development  of  the  virtual  work  equation,  a 
note  should  be  said  about  the  sign  convention  used 
here-  Compression  is  held  to  be  a  positive  stress  and 
tension  as  negative.  The  following  figures  illustrate 
the  sign  convention  used  here. 


> 


V 


I" 


> 


-Txz  <r 


A 


A 


<£ *ax 


Figure  (2-1)  Sign  convention  for  inplane  axial  and 

shear  stress  components. 
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■»  Y 


Figure   (2-2)    Application  of  sign  convention  for 

case  of    an  axial  load. 
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2-  3  Development  of  Virt  ua  1,  Work  Eguat ion._ 

Model  I  is  not  assumed  to  be  perfect  but  rather  to 
have  initial  imperfections.  These  initial  imperfections 
v  ,  w  and  B  and  the  additional  deformations  v,  w  and  B 
are  illustrated  in  figure  (£-3). 


Figure   (2-3)    Illustration  of   initial   imperfections 
v*,  w-*,  B*  and  additional  deformations  v,  w,  B. 


The  virtual  work  for  the  case  of  ari    inplane   axial 
load  becomes  (reference  1): 


W  -  Jj      5-   u(y,  z)  dy  dz 


A   x 

Thus  using  figure  (1-1),   it  can  be  shown  that  the 

strain   of   the   centroidal    axis    in   the    bent 

configuration,  e"   (membrane  strain),  is  (reference  1)2 

x 

2  2 

ft  -  -u   +  l/2(v  )   +  l/2(w  ) 

XXX  X 
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for   a   perfect   I   section.    With   imperfections,   we 

substitute: 
# 
v  +  v   for  v 

* 
w   *   w      for   w 

Thus: 

*  £  *   £ 

&      ■   -u      +    l/£(v   +v      )       +    l/2(w   +w      ) 

XX  XX  XX 

using  the  inextensibi lity  assumptions,  set: 

&   -  0 
x 
we  have: 

*  2  *  £ 

u      =    l/2(v  +v    )      +    l/£<w   +w    ) 
x  xx  xx 


f      C(v 

J  0\  i 


#£  *   2 

C(v   +v    )       +    (w   +w    )    ]    dx 
0  x      x  xx 


where   u  =   u(y, z) . 

From  geometrical  relationships  (see  figure  (2-3) ) 

v      *   zB  w      *   -yB 

♦  *  *  * 

v     *   zB  w      *  -yB 


(£-1) 


V         Sf     ZB  W         2!     -yB 

XX  XX 

*  *  ♦  * 

v      s*   zB  w      a*  ~yB 

XX  XX 

Substituting  and  integrating  we  get: 
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/■a      2   2  2        *  2      *2 

u(y,z)    =    1/2  J       (z    B      +   2z    B   B      +    z      B 

8  X  XX  X 


2   2  2         *2  2   *2 

+yB      +   2y   B   B         +yB      )dx 
x  xx  x 


r. 


2  2  *   2 

u(y,z)    =    1/2  (z      +   y      ) <B      +  B    )      dx 

x  x 


# 
Again  using  the  simple  choice  for  B  and  B  : 


Thus: 


B   - 

B   sin 
o 

(m  lr  x/a) 

* 

* 

B 

-   B   sin(m/?"x/a) 

o 

ra      2      2 

)    = 

1/2    J 

(z    +y    ) 

0 

*  2 

+    <(m^/a)B   cos(m/Vx/a)3      dx 
o 


/-a      2      2  *   2 

u(y, z)    ■    1/2    J      (z    +y    ) (B   +B    ) 
*M3  o      o 


2 
*    C<m  77Va)cos(m  Ifx/a)  1      dx 


2      2  *   2  2 

u(y,  z)    -    l/2(z    +y    )  (B   +B    )     (m^/a) 

o      o 


a 


*    Cx/2   +    (sin(2m^x/a)/(4rn/7'/a)3 

0 
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2   2      *  2       2 
u<y,z)  ■  1/4(2  +y  )  <B  +B  )  (m^/a)  a 

o   o 


The  work  equation  now  becomes: 

2      2  * 

:    +y    ) <B   +B    I 
Ax  o      o 


r  2      2  *   2  2 

W   =     \\       §5    Cl/4(z    +y    )  <B    +B    )     (m  7/Va)    a]    dy   dz 


For  the  current  situation  of  end  loading,  £   =  s" 

x    e 

*   2  2         rr       2      2 

W   =   §    C1/4(B   +B    )     (m  /77a)    all    j  J     (z   +y    )    dy   dz 
e  o      o  ^'A 


But    note   that: 


rr     2       2 

J  J       (z      +   y    )    dy   dz   =      I 
P.  P 


Thus: 


♦   2  2 

W  «   1/4  ft    (B   +B    )     (mtfVa)    a    I 
e      o     o  p 
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2. 4  Determination  of  the  Cri.ti.cai  Buckling  Stress^ 

Applying   the  calculus  of  variations  with   respect 

to   B    to  the  work  and  strain  energy  equations   yields 

o 
(reference  5) : 

*  2       2 
6   W  =  1/2  £  <B  +B  )  (m#Va>  a  I  &   B 
e   o  o  p   o 

2      2  2 

cfv  =  1/2  a  B  (m7>Va)  CE  <  I  s   +  C  )  (rn  7>/a) 

O  2         W 

+  gj:  i  b 

o 


Applying  the  principle  of  minimum  potential  energy: 

77"  =     £v  -  £w  =  0  , 

we   get : 

Solving    for   £    : 

e 

2  2 

CE(I    s      +   C    )  (m  /?7a>       +   GJ3  B 

z  w  o 


e  * 

I  B      +   B 

p  o           o 

(2-2) 

It   can   be   seen  by  inspection  that   the   lowest 

buckling   stress  occurs  for  one  wave,  m  =  1,   since   m 

must  be  an    integer.  Thus: 
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§5   -  §5 
e    CR 

As  it  can  easily  be  seen,   as  9   becomes  large  s" 

o  CR 

approaches  the  limit  of: 

£  2 

[Ed  s  +  C  )  (m  Tr /a)   +  GJ3 
z       w 

aj    « (2-3) 

CL 

I 
P 

Figure   (2-4)    illustrates   the  behavior  of  s"    with 

*  CR 

increasing  B   +  B  . 
o    o 


Further  simplification  and  rearranging  of  equations 
(2-2)  and  (2-3)  provides: 


.  B 
o 

8f     m    8{       

CR    CL         * 
B   +  B 
o    o 


SB+SB=SB 
CR   o  CR   o  CR   o 


* 
B 
o 

B      * (2-4) 


o 

S      /£        -    1 
CL      CR 
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B0  +  B0 


Figure    (2-4)       £        vs.     (B      +   B    ) 

CR  o  o 
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2. 5  Determination  of  t  he  First  Yield  Load._ 

To  determine  the  first  yield  load,  the  total  torque 
developed  in  the  torsional  buckling  process  must  be 
first  evaluated.  This  torque  will  consist  of  two  parts. 
The  first  term  consists  of  the  well  known  St.  Venant's 
Torque  (reference  1).  The  second  term  represents  the 
warping  torque  (reference  2). 
From  reference  ( 1 )  : 


T  ■  T    +  T 

SV     W 


where: 

T   ■  GJB 
SV      x     - 

and  from  reference  (2) : 

T  =  -  EI   d   w 

W  Z    C    XXX 

and  from  geometry; 


w   =  d  3 
X     ex 

resulting  in: 

2 
T  =  -  EI  d  B 

W  Z  C  XXX 

Thus  the  total  torque  is: 
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2 
T  =  GJB   -  EI  d  B  (2-5) 

X        2  C  XXX 

Figure  (2-3)  illustrates  the  above  equations. 


Again  a  simple  choice  of  B  is  made,  and  with  m  =  1: 

B  =  B  sinCtf-  x/a)  (2-6) 

o 

Equation  (2-5)  now  becomes: 

T   -   GJ(77Va)B   cos(/7-K/a) 
o 

2  3 

+   EI    d    (  7/Va)    B   cos<77-x/a) 
z   c  o 

2  2 

T   -    (TP-/a)B   cos(7>x/a)  CGJ   +   EI    d    (^/a>    1 

O  2  C 

Choosing  a  simple  sinusoidal  shape  for  T  allows: 

T  -  T  cos iV  x/a) 
o 
Thus: 

2      2 
T  =  (7?7a)B  CGJ  +  EI  d  </?Va)  1  (2-7) 

O  O  2  C 


To  eliminate  B  ,   equation  (2-4)  is  applied  to  equation 

o 
(2-7). 

* 
B 
o 
2      2 
T  -  CGJ  +  EI  d  (7^/a)  3("Va)  

O  2  C 

(§  /%-   -  1) 

CL   CR 

(2-3) 


37 


From  sectional  considerations,  an  evaluation  of  the 
maximum  stress  can  be  made. 


a 


X 


\ 


Figure  _(2-5>  Sectional  illustration. 

The   moment   developed  by  the  torque   in   terms   of 
deflections  is  (reference  2): 
M(x)  ■  -  EI  *  . 

2  XX 

So   the   maximum  compressive  stress   developed   in   the 
flange  due  to  the  torque  is: 


-EI  d  B   f 

2  C  XX  W 


(S  ) 

x  MflXT 


£1 


(S  )     -  -  Ed  B  f  /a 

x  MflXT        c  xx  m 


Application  of  the  above  definition  for  B  (equation 
2-4)  leaves: 


38 


Ed  f 

c  w  2 

(§!  )      - (/7"/a)  B  sin(/7-x/a) 

x  MftXT  o 

2 

A  choice  of  x  is  made  such  that  (s"  )      is  a  maximum: 

x  MftXT 

Ed  f 

c  w  2 

<£  > (  7^/a)  B  (2-9) 

x  MftXT  o 

2 


Equations  (2-7)  and  (2-9)  are  combined  to  eliminate  B  . 

o 

Ed  f  T  (  ^Va) 
c  w  o 

(§S  ) (2-13) 

x    MftXT  2  2 

2CGJ   +EI    d    (  tfVa)    3 

2  C 

Using  equation  (2-8)  leads  to:  m 

2        * 
Ed  f  <tfVa>        B 
c  w  o 

(£  )      = (2-11) 

x  MftXT 

2  Cs-   /£      13 

CL   CR 


Note  that   when  B   =0  that  this  equation  leads  to  art 

o 
indeterminate  form  of  5   since  s"   =  s* 

x        CL    CR 
The  total  maximum  stress  cart   now  be  described  by: 


(5  )       =  (£  )      +  P/fl 
x  TOTftL      x  MftXT 


where 
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P/A  ■  £ 


CR 


Thus  the  yield  stress  has  been  shown  to  be: 

§f   ■  (5  )      +  £ 
y      x  MflXT    CR 

Using  equation  (2-11)  in  equation  (2-12): 


(2-12) 


Ed  f  (77" /a) 

c  w 


B 


y 


£   /£    -  1 
CL   CR 


CR 


(2-13) 
2 
Solving   for  5    and   dividing  by  £   to   place   into 

CR  y 

nondimensional  form  yields: 


2  * 
S      Ed  f  (*fr7a)  B\   /S  \   /£ 

CR^      L-^+  ----- - --- W---U 

2%    /  \\  I  \\. 

(2-14) 
Choosing  the  smallest  root  and  simplifying  the  equation 
by  introducing  a  slenderness  ratio  L: 


L  - 


(2-15) 


\1 


CL 
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and  letting: 

2      * 
(  77  /a)  Ed  f  B 
c  w  o 

ft  -  (2-16) 

25 

y 

Then  : 

2 
5       (1  +  1/L   +  ft)  - 
CR 


2     2      2 

(1  +  1/L  +   ft)   -  4/L 


£  2 

y 

(2-17) 

ft   plot  of  s-   /§»  vs  L  for  various  values  of   ft   is 
CR   y 
shown   in  figure  (2-6).   Figure  (2-6)  was  compared  with 

various  design  criteria  from  reference  (10)  and  it   was 

found   that   equation   (2-17)   was   similar  to  design 

equations  of   reference   (10).   However,   the   reduced 

slenderness  ratio,   L,   of  reference  (10)  is  based  upon 

simple   Euler  buckling  and  the  geometric  parameter,   ft, 

is  of  course  different.   The  column  curves  of  reference 

(10)   were  quite  similar  in  appearance   and   range   of 

values  to  that  found  in  figure  (2-6). 
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Figure  (2-6)   Column  Curves 
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3.  MODEL  II  -  TORSIONAL  BUCKLING  UNDER  AXIftL  LOADS 

3.  1  Deveiogment  of  the  Strain  Energy  isuat ion._ 

Since  Model  II  is  the  very  simple  case  of  a  flat 
bar  stiffener,  thin  plate  theory  cari  be  easily  utilized 
in  the  development  of  the  strain  energy  equation.  Much 
of  this  section  has  been  patterned  after  a  similar 
analysis  presented  in  reference  (3). 

For  a  web  plate  of  uniform  thickness  and  a  length  a 
the  strain  energy  is  repensented  by  (reference  5) : 


•3D  f*f 


<x)  2 

D   )   /      C(v    +  v   ) 

W  ^  O     'O  XX      ZZ 

2 
2(1-  N) Cv   v   -  v   ]]  dx  dz 
XX  zz     xz 


f 


a   2 
CB   dx 
o 


The   second   term  represents  the   rotation   of   the 

supporting   plate  structure  modeled  here  as  an  elastic 

spring.  D   is  the  flexural  rigidity  of  the  web  plate. 

w 
From  geometrical  considerations: 

v  ■  zB 

making  a  simple  choice  for  B; 

B  ■  B  sin  (m  7?-x/a> 
o 
substituting   into  the  expression  for  the  strain  energy 

and  integrating  over  a  constant  depth  d  results  in: 
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2   2  2  2 

V   -   D      ad    (rnTTVa)    B       C    d     (m/^/a)       +   6(1    -V)    ]/i2 

w  o 


+   CB   a/4 
o 


(3-1) 
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3.  2  Deyel^ogment  of  the  Vi.rt  ual^  Work  isuat  i.gn._ 

As   for   Model   I,   Model  II  is  not  assumed   to   be 

perfect  but  rather  to  have  initial  imperfections.  These 

*     *     * 
initial   imperfections  v  ,   w  ,   B   and  the   additional 

deformations   v,   w,   and  B  are  illustrated   in   figure 

(3-1). 


Figure  (3-1)   Illustration  of  the  initial  imperfections 
v*,  w*,  B*  and  additional  deformations  v,  w,  B. 


For   the   case  of  an  inplane  axial  load,   the   work 
equation  is  : 


W  ■  jj   S  u(y,  2)  dy  dz 


From  figure  (3-1),  it  can  be  shown  that  : 

2  2 

&      -   -   u      +    l/2(v    )       +    l/2(w    ) 

XXX  X 
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for   a  perfect  web   section.   With   imperfections,   and 

using  the  inextensibi 1 ity  assumption  with  regard  to  the 

length  of  the  stiffener  setting  e*   =  0,  the  result  is  : 

x 

#2  #2 

u      =    1/2    (v      +   v    )       +    l/2(w      +   w    ) 
x  xx  xx 


a  *   2  *   2 

C(v      +   v    )       +    l/2(w      +   w    )    1    dx 

O  X  X  XX 

(3-2) 


where   u  =   u(y, z) 

As   shown   in   figure    (3-1),    the    geometrical 

relationships   (equations   2-1)   are  used   to   simplify 

equat  ion  (3-2) : 

fa        2    2       * 
u(y,z)  =  1/2   J    (z   +  y  ) (B  +  B  )  dx 
'  o 

* 
A  simple  choice  is  made  for  B  and  B  . 

B  =  B  sin  (m>x/a) 
o 

*    * 
B   =  B  sin  (m/7"x/a) 
o 

Thus: 

2    2        *        2 
u(y,z)  =  l/4(z   +  y  )  (B   +  B  )(mJ^/a)  a 

o    o 

The  work  equation  now  becomes: 


w  =  ff»  \ 


2    2        *        2 
Cl/4(z   +  y  ) (B   +  B  >(mff-/a)  a3  dy  dz 

o    o 
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*  2       2 
W  =   1/4  3  (B   +  B  )  (m//Va)  a  I 
x   o    o  p 


3. 3  Determination  of  the  Crit igai  Buckling  St ^ess^ 

Applying  the  calculus  of  variations  with  respect  to 

B   to  the  work  and  strain  energy  equations  yields: 
o 

*   2  2         f 

a  W   =      1/2   £    (B      +   B    )     (m  /7"/a)    al    d  B 
x      o  o  p        o 

r  2  2  2 

6  V   ■   -C    D      ad    (m^Ya)    B      C   d    (m  ^/a) 


w 


2  r 

+   6(1    -V)    3/6      +   CB    a/2    >dB 


Applying  the  principle  of  minimum  potential  energy: 

77  s  i"v  -  Sw  =  8 

<^W  =  §  V 

Solving  for  s"  ■  3   : 
e    CR 

2      2  2 

D   d(m77Va>     Cd    <m  7^/a)       +   6(l-V)3/3   +   C  B 

w  o 

CR  2  * 

(m  #7a)  I  B      +   B 

p  o  o 


(3-3) 
For   the  case  of  zero  restraint   against   rotation; 

C  ■  0: 
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2  2       2 

D  d(m^/a)  Cd  (m>/a)   +  6(1-V)3        B 
w  o 

£     a 

CR  * 

3  1  B   +  B 

p  o     o 

(3-4) 

It  is  easily  seen  that  the  lowest  £    for  this  case 

CR 
results  from  the  simplest  of  mode  shapes,  rn=l. 

When  the  rotational  restraint  is  not  equal  to  zero; 

the   lowest   value  of  5*   depends  upon   the   degree   of 

CR 
rotational  restraint: 

3  1/4 
m  ■  <a/7f ) (3C/D  d  )  (3-5) 

w 

With  the  use  of  equation  (3-5),   m  may  be  evaluated 

an,d   rounded  up  or  down  to  the  nearest  integer  for   the 

correct  wave  shape. 

It   is   easily   seen  that  as  B   becomes   large   £ 

o  CR 

approaches  a  limit.  Thus  with  (s"   )       -   s"   : 

CR  LIMIT     CL 


2  2       2 
D  ddn^/a)  Cd  (m7?-/a)   +  6(l-V)3/3  +  C 
w 

£    « 

CL  2 

(m  "fr/a)     I 


(3-6) 

Figure   (3-2)   illustrates   the  behavior   of  %"    with 

*  CR 

increasing  B   +  B  . 
o    o 
Further   simplification   and    rearranging     of 
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equations  (3-3)  and  (3-6)  provides; 

* 


B 


B   ■ 
o 


5   /S   -  1 
CL   CR 


(3-7) 


'CR 


A 


CL 


B, 


B0  +  B0 


Figure  (3-2)   Behavior  of  S   with  B   +  B 

CR       o    o 


The  first  yield  load  will  not  be  determined  for 
this  case  under  this  method  of  analysis  due  to  the 
relative  complexity  of  equation  13-3)  and  due  to  the 
fact  that  the  Perry-Robertson  method  will  not  lend 
itself  to  further  simplification  along  this  train  of 
analysis.  Since  a  simplified  formula  for  design 
purposes  is  the  goal  of  this  thesis,  the  next  chapter 
will  use  a  different  tactic  to  solve  the  problem. 
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4.  MODEL  II  -  GENERALIZED  ANALYSIS 

4.  1  Method  of  Approach  and  Strain  Energy  H9y§£i°.Gs. 

The  analysis  thus  far  for  the  case  of  a  flat  bar 
stiffener  under  end  loading  has  proceeded  along  the 
lines  of  thin  plate  theory  up  to  the  development  of  art 
expression  for  the  critical  buckling  stress. 

This  chapter  uses  beam  theory  to  develop  a  simpler 
formulation  in  a  manner  analagous  to  that  performed  on 
Model  I  in  chapter  £.  In  fact,  this  analysis  will 
merely  be  a  generalization  of  the  results  for  Model  I 
applied  to  Model  II.  It  must  be  remembered  here  that 
Model  I  consisted  of  two  unequal  flanges  with  the  top 
flange  deforming  or  rotating  about  the  bottom  flange. 
At  no  time  in  the  derivation  of  the  formula  for  the 
critical  stress  of  the  tee  stiffener  is  there  an 
assumption  made  about  the  position  of  the  top  flange 
relative  to  the  bottom  flange.  Thus  the  same  analysis 
with  minor  changes  will  be  performed  for  the  flat  bar 
stiffener.  However,  in  the  case  of  the  flat  bar 
stiffener  of  Model  II,  there  is  a  restraining  effect 
provided  by  the  juncture  of  the  flat  bar  with  the 
plate.  This  in  effect  can  be  modeled  as  a  spring  with  a 
rotational  constant  C.  Thus  the  strain  energy  equation 
would   have  four  terms  (instead  of  the  three  terms   for 
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Model  I).  These  terms  represent  (in  order  of 
appearance)  sideways  bending;  longitudinal  warping; 
torsion;  and  rotation  of  the  supporting  plate  structure 
modeled  as  an  elastic  spring.  Thus  the  strain  energy 
equation  is: 


■»f 


"a      £        2       £     2 
J    CEI  v    +  EC  B    +  GJB   +  CB  1    dx 

O        2  XX        W  XX         X 


Again  choose: 

B  ■  B  sin<m  rJ"x/a) 
o 

Substituting  and  integrating  results  in: 

2        2  4  2 

V  =  1/4  a  B   C  E(I  s   +  C  )  <m  ^Va>   +  GJ(m7/7a)   +  C  1 
o      z      w 

Applying   the   calculus  of  variations  with   respect   to 

B   yields: 
o 

2  4 

S  V  m    1/2  a  B   C  E(I  s   +  C  )  (ra  i'/a) 
o      z      w 

B 

o 


+  GJ(m  "Tf/a)   +  C  3  o> 


4. 2  The  Work  Eguat ion.. 

The  assumptions  used  in  deriving  the  work   equation 

for  Model   I  are  the  same  as  for  Model  II.   Thus   from 

chapter  2  we  have: 

♦  2       2 
W  -  1/4  S   <B   +  B  )  (m  v?7a)   a  I 
e    o    o  p 
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and  r 

S   W  =  1/2  S   (B   +  B*>  <m  7?-/a)2  a  I   ^  B 


(4-2) 


4.  3  Det erminat  ign  of  the  Critical,  Buckiing  St ress 

Using   equations  (4-1)  and  (4-2)  and   applying   the 

principle   of  minimum  potential  energy  and  solving   for 

the  stress  yields  (with  s"   =  £   )  : 

e    CR 

2  4  2 

E(I  s   +  C  )(m7>Va>   +  GJ(m7T/a)   +  C   B 
z       w  o 

CR  2  * 

(m^-/a)   I  B   +  B 

p  o    o 

(4-3) 

Lets 

2  4  2 

E(I  s   +  C  )(m#Va)   +  GJ(m"ft*/a)   +  C 

Z         W 

5    = (4-4) 

CL  2 

(m  #7a)   I 

P 
Using   (4-4)   in   equation  (4-3)   and   solving   for   B 

o 
yields: 

* 

B 
o 

B (4-5) 

o 

£  /%        -    1 
CL  CR 

Comparison   of   equation  (4-3)  with  equation   (2-1) 
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reveals  the  only  difference  is  that  equation  (4-3)  has 
an  additional  term  in  the  numerator  (specifically  the 
spring  constant  C). 

Since  it  is  very  difficult  to  quantify  the  spring 
constant  C  ar\d  (hopefully)  since  its  value  is  small  in 
comparison  to  other  terms  in  expression  (4-3),  C  will 
be  set  to  zero  value  and  ignored  for  the  remainder  of 
this  analysis. 

Thus  for  the  case  of  C  ■  0,  equation  (4-3)  reduces 
to  equation  (2-2). 


4. 4  Determination  of  the  First  Yiel,d  Load.. 

The  determination  of  first  yield  load  for  the 
generalized  case  proceeds  in  a  similar  fashion  as  for 
the  case  of  Model  I. 

Again  the  total  torque  consists  of  two  parts 
(reference  1),  St.  Venant's  torque  and  a  warping 
torque.  However,  in  this  case  the  warping  torque  is 
(reference  2) : 


T   ■  -EI  d  w 

W         Z     XXX 


and 


w   =  dB 
x      x 


resulting  in: 
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2 
T  -  -EI  d  B 

W  Z     XXX 

Thus  the  total  torque  is: 

2 
T  ■  GJB   -  EI  d  B 

X        Z     XXX 

Choosing  B  and  T: 

B   -   B   sin(m7fx/a) 
o 

T   3   T   cos(m7?"x/a) 
o 
Leaves: 

2  2 

T      =    <7?Va>CGJ   +   EI    d    (  TTVa)    3    B  (4-6) 

o  z  o 

Using  equation  (4-5)  to  eliminate  B   leaves: 

o 

* 
B 
2      2         o 

T     ■    0?/a)[GJ   +   EI    d    (  7?7a)    1    

o  z  « 

(§5       /£  -     1) 

CL      CR 

(4-7) 

Again  sectional   considerations   can   be   used   to 

evaluate   the   maximum   stress.    The   maximum   stress 

developed  in  the  stiffener  due  to  the  torque  is  now: 

2 
(5  )      -  -Ed  B 
x  MAXT        xx 

Application   of   the   above  definition  of  B   and   a 

choice  of  x  such  that  (£  )      is  a  maximum  leaves: 

x  MAXT 

2      2 
(s-  )      =  Ed  (  /?7a)  B 
x  MAXT  o 
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Using  equation  (4-6)  to  eliminate  B  : 

o 
2 
Ed  T  (  T/a) 
o 
<§  )      s (4-8) 

x  MflXT  2      2 

CGJ  +  EI  d  i^/a)    1 

2 

Using  equation  (4-7)  in  equation  (4-8)  leads  to: 

* 
B 
2      2       o 

(gj  )      =  Ed  <y?7a)   (4-9) 

x  MflXT 

<£   /S    -  1) 
CL   CR 

Again,   from   sectional  considerations,   the   total 

maximum  stress  can  be  described  by: 

S  =  (5  )       +  £  (4-10) 

y      x  MflXT     CR 

Using  equation  (4-9)  in  equation  (4-10),  solving  for 

2 
5"    and   dividing  by  s"   to  place   into   nondimensional 

CR  y 

form  yields: 


2      2  * 
Ed  (~P"/a)  B  \  /U 


o 


CR1 


y        y 

(4-11) 
Choosing   the   smallest   root  and   simplifying   the 
equation   by  the  use  of  the  reduced   slenderness   ratio 
(equation  2-15)  but  redefining  parameter  A  as: 
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2   2* 
(  ^Va)  Ed  B 


A  = 


(4-12) 


Results  in  (equation  2-17): 
3 


2 


(1  +  1/L   +  ft)  - 


CR 


N 


2     2      2 
(1  +  1/L   +  ft)   -  4/L 


£ 


(4-13) 
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5.  COMPARISONS  WITH  OTHER  STUDIES. 

5. 1   Comgar ison   between   Th in  Plate   and   Beam   Ihegry, 
Results 

In  chapter   3,   the  critical  buckling  stress  of   a 

flat  bar  stiffener  was  derived  through  the  use  of   thin 

plate  theory: 

2   2       2 
D  d(m77Va)  Cd  (m^/a)   +  6(l-V)3/3  +  C    B 
w  o 

CR  2  * 

(m  ft/a)       I  B   +  B 

p  o    o 

(3-3) 

In  chapter   4,   the  critical  buckling  stress  of   a 

flat  bar  stiffener  was  found  using  one  dimensional  beam 

theory: 

2  4  2 

E(I    s      +   C    MrnT/Va)      +   GJ(m"Va)       +   C  B 

z  w  o 

£  . 

CR  2  * 

(m^/a)         I  B      +   B 

p  o  o 

(4-3) 
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It  should  be  noted  here  that  when: 

2      2  2 

D   d<m/7*/a>     Cd    <m  7?7a>       +   6(l-V)]/3      - 
w 

2  4  2 

E(I    s      +   C    )  (m  TT/a)       +   GJ(m7?Va) 

z  w 

the  agreement  between  the  two  methods  would  be  exact. 

For  the  elastic  range  of  experimentation  and  with: 
3 
Et 
w 

D   = G  *  E/2(l+V)  =  E/2. 6 

w 

12U-V) 


3  3 

I   -  <d  t  )/12  J  =  (d  t  )/3 

2        W  W  WW 

d  =  d  V  =  0.3 

w 

The  above  relationship  reduces  to: 


2       2  2       2 

d  (m7/Va)  1         d  (m7TVa) 

w  w 

+     9     


10.92  2.6  12  7.8 

(Thin  Plate  Theory)  (Beam  Theory) 

From  a  cursory  examination  it  would  appear  that  the 
Thin  Plate  Theory  results  (equation  3-3)  will  probably 
give  a  larger  value  for  the  buckling  stress  for  the 
normal  range  of  stiffener  dimensions  used  in  ships.  To 
test   this   hypothesis,   the  dimensions  of  a   flat   bar 
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stiffener   used  in  experiment  #21  of  reference  (6)  will 

be  used. 

where: 

d   ■  87. 4  cm 
w 

a  =  730  cm 

m  -  1 

For  the  thin  plate  theory: 

£       2 
d  (m7T/a)  1 

w 

+    =   0.398 

10.92  2.6 

For   the    beam   theory: 

2  2 

d    (m  TT/a)  1 

w 
+        =      0.  14 

12  7.8 

Thus,  within  the  elastic  range,  the  thin  plate 
theory  results  would  indeed  provide  higher  buckling 
stress.  The  beam  thoery  results  may  thus  be  too 
conservative.  However,  section  5.3  shows  good  agreement 
between  equation  (4-11)  (first  yield  load)  and 
experimental  yield  results. 


5. 2  Comparisons  wi,t h  Published  Formulae;, 

Reference  (3)  performed  a  similar  analysis  to  that 
performed  in  this  thesis  but  without  the  inclusion  of 
initial  geometric  imperfections. 

For   the   case  of  a  Tee  stiffener  subjected   to   an 

axial  load,  equation  21  of  reference  (3)  allows  (in  the 

notation  of  this  thesis): 

2  2  2 

GJ  +  Ed  s   +  C  Xm7?Va>   +  C(a/m<f) 

z      w 

CR 

I 
P 

Equation  (2-2)  of  this  thesis,   for  a  Tee  stiffener 

in  an   analysis  which  ignores  the  contribution  of  the 

web   and  includes  the  contributions  of  initial  defects, 

is: 

2  2 

[Ed  s   +  C  )  (m  777a)   +  GJ3        B 
z      w  o 

S 

CR  ♦ 

I  B   +  B 

p  o    o 

Thus   equation  (2-2),   for  the  case  of  an  initially 

perfect   Tee   stiffener  reduces   to   equation   21    of 

reference   (3)  when  the  rotational  restraint  is  held  to 
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be   zero.   However,   it   must  be  noted  that   since   the 

physical   models  of  the  two  derivations   were   slightly 

different  -  namely  that  the  model  used  in  this  analysis 

consisted    of   two   unequal   flanges   -  some   of    the 

parameters   of  the  equations  will  have  different  values 

(ie.  I  ,  J  etc. ) . 
z 
For  the  case  of  a  flat  bar  stiffener,   equation   £2 

of  reference    (3)  has: 

2  2  2 

D  dCCm^/a)  d   +  6<1-V)3       C(a/m/T) 
w 

*. +    

CR 

31  I 

P  P 

From   chapter    4,     equation    (4-3): 

2  4  2 

E(I    s      +   C    )  (m  T/Va)       +   GJ(mfiVa)       +   C  B 

z  w  o 

5-  = 

CR  2  * 

I     (mTTVa)  B      +   B 

p  o  o 

The  following  is  a  comparison  for  a  Tee  stiffener 
and  a  flat  bar  stiffener  which  was  presented  in 
reference  (3).  Figure  (5-1)  from  reference  (3)  shows 
the  geometric  dimensions  of  the  Tee  and  flat  bar 
stiffeners.  Table  (5-1)  shows  how  equations  21  and  22 
of  reference  (3)  and  equations  (2-2)  and  (4-3)  of  this 
thesis  compare  with  a  finite  element  analysis  performed 
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using   a   computer  program  developed  and  documented   at 
the  University  of  California  (reference  12). 


4 
Flat  Bar         (5   /E)  x  10          C  =  0 

CR 

FEM 

EQN<22) 

EQNC4-3) 

m  -  1 

6.44 

6.42 

6.26 

m  -  2 

7.58 

7.51 

7.03 

4 
Tee              <£   /E)  x  10          C  -  0 

CR 

FEM 

EQN(21) 

EQNC2-2) 

m  ■  1 

26.4 

26.3 

30.69 

m  =  2 

65.9 

72.6 

89.4 

Table  (5-1)  Comparison  with  Published  Formulae 
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t^  -  0.635  cm 


d  ■  d-  16  cm 


t  -  0S4  cm 


b#  -  30  cm 


FLAT  BAR  STIFFENER  -  a  -  100  em 


1W  -  7.90  cm 


^  •  1.42  cm 

tw  -  0.722  cm » 

N 


t  -  0.80  cm 


d  •  15.75  cm 


b  •  S.4"»  cm 

i 

_I 


14.33  cm 


f 
I 

b§  -  30.5  cm 


TEE  STIFFENER  -  a  -  160  cm 


Figure  (5-1)   Stiffener  Geometries  for 

Comparative  Solutions 
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The   U.S.   Navy   has   in  the   past   used   similar 

formula   as  those  presented  here  for  the   determination 

of   critical   buckling   stresses.    For   the   critical 

buckling  stress  of  a  Tee  stiffener,   in  the  notation  of 

this  thesis  from  reference  13  we  have: 

2 
GJ  +  (rntt/a)  EC 

w 

5    = 

CR 

I 
P 

where  here: 

3  2         3  3 
C   -  f  d  t  /12  +  d  t 
w     w    f         w 

3        2        3 
I   =dt/3+dft   +ft 
p       w         w  w     w  f 

This  is  very  similar  to  those  results  determined  by 
reference  (3),  with  the  exception  that  the  above 
formula   does  not  account  for  the  rotational   restraint 

provided  by  the  junction  of  the  web  and  plate.   Also  I 

P 
is  much   simplified  here  and  C  ,   the   torsion-bending 

w  2 

constant,   equates   to   the   quantity  <I  s   +  C  )   of 

z  w 

reference  (3)  -  but  also  somewhat  simplified. 
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5.  3  Comparisons  with  £><fier  i.ment a.1  Resu  lts._ 

Reference  (6)  provides  insight  into  the  torsional 
buckling  process  through  the  experimental  results 
published  there.  This  report  describes  model  tests  on 
the  collapse  and  post  failure  strength  of  Tee  struts 
simulating  flat  bar  stiffeners  in  a  stiffened  plate. 
The  T  struts  were  subjected  to  an  axial  load  up  to  and 
beyond  failure.  The  geometrical  imperfections  of  the 
models  were  also  recorded.  The  researchers  of  reference 
(6)  (of  Det  norske  Veritas)  are  to  be  lauded  for  the 
thoroughness  and  completeness  of  their  report, 
particularly  in  the  treatment  of  the  experimental  data. 

Table  (5-2)  contains  the  model  data  pertaining  to 
this  analysis  for  ten  experiments  of  reference  (£). 
Figures  (5-2  to  5-11)  show  the  load-  deflection  curves 
for  these  ten  experiments.  The  solid  line  represents 
the  experimental  data.  Point  ft  of  each  of  these  graphs 
represents  the  first  yield  load  point  as  determined  by 
eq  uat  i  on  (4-11). 

Notice  that  there  is  very  good  correlation  between 
the  experimental  first  yield  load  point  and  the  results 
of  equation  (4-11). 
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Mode 

L    a 

h 

t 
w 

# 
B 
o 

y 

CR 

u 

21 

703 

87.  4 

4.0 

1.67E-02 

256 

233 

230 

22 

703 

87.3 

4.0 

5. 6E-03 

271 

263 

246 

23 

703 

87.7 

4.0 

1. 15E-02 

273 

257 

234 

24 

703 

87.7 

4.0 

6. 0E-03 

271 

262 

240 

25 

703 

87.7 

4.0 

8. 0E-03 

272 

260 

230 

26 

703 

87.7 

4.0 

8. 0E-03 

256 

244 

243 

27 

703 

87.4 

4.0 

6. 0E-03 

266 

257 

272 

28 

703 

88.0 

4.0 

7. 99E-03 

256 

244 

241 

29 

703 

87.8 

4.0 

4. 96E-03 

259 

252 

218 

3)3 

703 

88.  1 

4.0 

1.26E-02 

271 

253 

230 

Not at  ion: 

a  -  length  of  model  (cm) 

h  -  height  of  stiffener  (cm) 

t   -  thickness  of  stiffener  (cm) 

M 


B   -  angle  of  rotation  (radians) 
o 


S  -  yield  stress  (N/mm  ) 

y 

s*   -  first  yield  load  stress  (N/mm  ) 

CR 

2 

S  -  ultimate  stress  (experimental  results) (N/mm  ) 
u 

Table  (5-2)  Data  of  Model  Tests 
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Figure  (5-2)   Model  Test  #21  of  Reference  (6) 
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Figure  (5-3)   Model  Test  #22  of  Reference  (6) 
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Figure  (5-4)   Model  Test  #23  of  Reference  (6) 
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Figure  (5-5)   Model  Test  #24  of  Reference  (6) 
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Figure  (5-6)   Model  Test  #25  of  Reference  <6> 
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Figure  (5-7)   Model  Test  #26  of  Reference  <6> 


72 


C   02 

J  .  00 


2.70 


2.40 


2.10 


1.80 


n    1.50 

ii 

II     1.20 

i 

M 

'•I     0.90 


0.60 


1          1 

i 

1 

h 

i 

\ 

\ 

1 

y 

1 

/ 

/ 

^~_ 

/ 

/ 

0.30 


0.00 

0.00       0.70       1.40       2.10       2.80       3.50       4.20       4.90       S.bO       o.jO       7  .few 


STP0KE<nf1) 


t    W 


MODEL    27 


Figure  (5-8)   Model  Test  #27  of  Reference  (6) 
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Figure  <3-9>   Model  Test  #28  of  Reference  (6) 
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Figure  (5-10)   Model  Test  #29  of  Reference  (S) 


75 


K    &2 

J.  00 


2.70 

2.4* 

ft 

5     1  -80 


ii    1.5C' 

11 

II    1.20 

r 

M 

'••    0.90 


0.60 


0.30 


0 .  00 


\     1 

1 

1 

A« 

EXPERIMENT 

A 

n 

. 

/ 

\ 

/ 

l 

0.00       0.70       \.4Q       2.10       2.80       3.50      «4.20       -t.yo      5.60       6.3y       V.Om 

STROKE  (fin)  L    H" 

MODEL     30 


Figure  (5-11)   Model  Test  #30  of  Reference  (6) 
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6.  CONCLUSIONS 

This  thesis  utilized  a  design  approach  to  tackle 
the  problem  of  torsional  buckling  of  stiffeners  with  a 
good  deal  of  success.  Integrated  into  all  phases  of 
this  analysis  was  a  consideration  of  the  geometric 
imperfections  of  the  stiffener  -  both  Tee  and  flat  bar. 
These  structural  members  were  analyzed  as  simplified 
sections  and  the  Perry-Robertson  approach  used  in  other 
areas  of  engineering  was  successfully  applied  here. 

Chapter  5  showed  that  the  derived  results  of  this 
thesis  reduced  to  currently  published  formulae  for  the 
case  of  the  perfect  stiffened.  The  contribution  of  this 
thesis  is  that  the  initial  imperfections  were  taken 
into  consideration,  using  a  simple  design  approach,  in 
the  derivation  of  the  results  presented  here. 

Also  shown  in  chapter  5  was  a  generally  good 
correlation  between  results  of  this  thesis  and  axial 
load  tests  on  flat  bar  stiffeners.  Comparison  of  the 
Tee  stiffener  formula  with  experimental  data,  which 
measures  and  reports  initial  imperfections,  needs  to  be 
performed.  Also,  comparisons  should  be  made  with  other 
tests  on  flat  bar  stiffeners  in  more  dimensional  ranges 
to  determine  the  range  of  validity  for  the  results 
derived  here. 
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Further  work  needs  to  be  done  in  this  area. 
Specifically,  determination  of  the  critical  buckling 
stress  and  first  yield  load  due  to  lateral  loadings. 
Also,  a  simplified  method  of  handling  the  mode 
interactions  of  combined  lateral  and  axial  loadings  is 
needed.  In  addition,  it  is  an  asset  to  understand  from 
the  work  reported  in  references  (6)  and  (15)  that  the 
falling  path  of  the  load-deflection  curve  (figure  6-1) 
is  represented  by  the  tripping  mechanism  which  is  not 
covered  in  this  analysis. 
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Figure  (6-1)   Experimental  Load  Deflection  Curve 
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